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ABSTRACT 

Wc report nine new transit epochs of the extrasolar planet WASP-5b, observed in the Bessell 
I band with SOAR at the Cerro Pachon Observatory and with the SMARTS 1-m Telescope at 
CTIC0 , between August 2008 and October 2009. The new transits have been combined with 
all previously published transit data for this planet to provide a new Transit Timing Variations 
(TTVs) analysis of its orbit. We find no evidence of TTVs RMS variations larger than 1 min 
over a 3 year time span. This result discards the presence of planets more massive than about 
5 , 1 M® and 2 M ffi around the 1:2, 5:3 and 2:1 orbital resonances. These new detection 
limits exceed by ^ 5 — 30 times the limits imposed by current radial velocity observations in the 
Mean Motion Resonances of this system. Our search for the variation of other parameters, such 
as orbital inclination and transit depth also yields negative results over the total time span of the 
transit observations. This result supports formation theories that predict a paucity of planetary 
companions to Hot Jupiters. 

Subject headings: exoplanets: general — transiting exoplanets: individual(WASP-5b) 



1. Introduction 

Once the method of Transit Timing Variations 
(TTVs) was theoretical proposed as of great po- 
tential to dete ct additional exoplanets in tran- 
siting systems dMiralda-Escudej 2002 : Agol et al 
20051 : iHolman fc Murray! I2005T ). and even ex- 



omoo ns ( Sa rtoretti fc Schneider 119991 : iKippinsJ 
20091) . several observational groups have started 
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to monitor the majority of the known transit- 
ing planets. This monitoring aims at detecting 
changes in the predicted mid-time of the transits 
to infer the presence of additional planets in the 
system not detected previously by e.g. radial ve- 
locities. Those data have also been used to detect 
variations in other transit parameters (e.g. tran- 
sit depth and duration), that can be attributed 
to perturbations produ ced by unseen companions 
(|Miralda-Escudell2002t) . 

In addition to having the potential of finding 
planets in the Earth-like or smaller mass regime, 
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the detection (or non-detection) of companions of 
transiting Hot Jupiters through TTVs also can im- 
prove constraints on planet formation models (e.g 
Triaud et ai]|20ld: iNaoz et alj|201lt iMiguel et al 



Gillon et al. (2009) was base d on the divergence 



201 LL and references therein) and help discrimi- 
nate between the different mechanisms proposed. 
In this way, TTVs become a powerful tool for the 
detection and study of multi-planet system ar chi- 
tectures (|Latham et al.ll201ll iFord et al.ll201lh . 

The first u nquestionable evidenc e of TTVs was 
announced by iHolman et al.l ( 2010t ) in the double 
Saturn- like transiting planetary system Kepler-9, 
where the central times of transit vary with ampli- 
tudes of 4 and 39 minutes in timescales of about 
19 and 40 days, respectively. Another extraordi- 
nary confirmation of the T TVs effect came wit h 
the discovery of Kepler- 11 (jLissauer et al.ll2011l) . 
a system with six transiting planets which shows 
TTVs of amplitudes as large as tens of minutes 
produced by the gravitational perturbations be- 
tween the planets. An additional remarkable re- 
sult of this later work has been the use of dynam- 
ical orbital fits to determine directly the masses 
of the transiting planets, dismissing the need of 
radial velocities. 

In 2008, the WASP-South survey reported their 
second detection of an exoplanet, WASP-5b, tran- 
siting a relatively bri ght star (V = 12.3) in the 
Southern Hemisphere ( Anderson et al.ll2008 , here- 
after A08). This discovery paper, based on WASP 
photometry and two additional transit epochs plus 
radial velocities measurements, announced a Hot- 
Jupiter planet with a mass of Mp = 1.58 



+0.13 
^0.08 



Mi 



and a density of p p = 1.22_ ; 24 pj, orbiting a G4V 
star with a period of P = 1.62 days. 

Gillon et all (|2009h did a reanalysis of the A08 



data to produce the first timing study of WASP-5b 
and arrived to the conclusion of potential period 
variations, based ona~ 2-minutc shift in the tim- 
ing residuals of the most precise points. 

l|2009h . 



Southworth et al 



hereafter S09, ob- 
served two new transits, for which they achieved 
very high photometric precision by defocusing 
the images at the 1.54-m Danish Telescope at 
La Silla Observatory, but at the expense of pro- 
ducing only 3-minute cadence. They refined the 
linear ephemeris of the system and concluded the 
high deviation of the timing residuals with re- 
spect to that straight line (Xred = 5-7) found by 



Smith et all (|2009l ) 



of only one point out of six. 
searched for signatures of additional planets in the 
residuals of WASP light curves after removing the 
transits of WASP-5b, and found no evidence of a 
transiting companion down to Saturn-size planets 
within periods of up to 20 days. 

Other recent works have determined and refined 
several physical parameters of the system. For ex- 
ample, iTriaud et all (j2010l) determined the angle 
between the orbital plane of WASP-5b and the 
spin axis direction of its host star to be consistent 
with zero (A = 12°ig°). This conclusion has been 



confirmed by the reanalysis of lFukui et al.l (|201ll ) , 
hereafter Fll, who obtain A = 7.2° ± 9.5. 

Fll additionally searched for TTVs of WASP- 
5b using seven new transit epochs, combined with 
all previously available observations. They find a 
RMS of about 68 seconds in their timing resid- 
uals despite of having an average of 41 seconds 
uncertainty per epoch, and proposed that such a 
large deviation from a linear fit (x 2 = 32.2 for 9 
degrees of freedom) can be explained by an orbital 
perturbcr. Using dynamical simulations Fll con- 
strained the masses of this hypothetical perturber 
to 2 M® in the 1:2 and 2:1 mean motion reso- 
nances (MMRs) and set a mass of 43 M® for a 
potential Trojan body. 

Dragomir et al.l ( 201lh . hereafter Dll, reported 



two new transits of WASP-5b with data of the 1-m 
telescope at Cerro Tololo Inter- American Observa- 
tory. 

In this work we present nine additional transits 
of WASP-5b, observed between August 2008 and 
October 2009, and perform a new homogeneous 
timing analysis of all available epochs to further 
confirm or rule out the TTV signals previously 
proposed for this system. 

In section [2] we describe the new observations 
and the data reduction. Section[3]details the mod- 
eling of the light curves and in section|4]we present 
the timing analysis. In Section [5] we discuss the 
mass limits for a unseen perturber. Finally, we 
present our conclusions in section [6] 

2. Observations and data reduction 

In 2008 we started the Transit Monitoring in 
the South Project, which is a monitoring campaign 
of transiting planets observable from the Southern 
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Hemisphere (jHover et al.ll201ll ). following the ap- 
proach of using high-cadence observations and the 
same instruments and setups to try to minimize 
systematics and reduce uncertainties in the mid- 
transit times, as well as other transit parameters. 
For the TraMoS project we have already observed 
more than 60 transits of over 20 exoplanets. 

As part of TraMoS we observed a total of 
nine transits of WASP-5b, between August 2008 
and October 2009^, with the Y4KCam on the 
SMARTS 1-m Telecope at Cerro Tololo Inter- 
American Observatory (CTIO) and with the 
SOAR Optical Imager (SOI) at the 4.2-meter 
Southern Astrophysical Research (SOAR) tele- 
scope in Cerro Pachon. 

Y4KCam is a 4064 x 4064 CCD camera with a 
Field of View (FoV) of 20 x 20 squared arcmin- 
utes and a pixel scale of 0.289 arcsec pixel -1 . The 
standard readout time of the camera is 46 sec, 
which we reduce to ~ 16 sec by binning 2x2. The 
SOI detector is composed of two E2V mosaics of 
4096 x 4096 pixels with a scale of 0.077 arcsec 
pixel -1 , giving a FoV of 5.2 x 5.2 squared arcmin- 
utes. The instrument has a 20.6 sec standard read- 
out, which becomes only ~ 11 sec after binning 
2x2. 

All nine transits were observed using a Bessell 
I filter (A cff = 8665 A and FWHM=3914 A) to 
reduce limb darkening effects in our light curves. 
Six of the transits were fully covered in phase. 
A fraction of the ingress of the 2008-11-03 tran- 
sit was not observed because a telescope system 
crash as illustrated in Figure [T] Nevertheless, this 
transit was treated as a complete transit. Two 
other transits, 2009-08-05 and 2009-10-21, were 
only partially covered with data between phases 
-0.034 < (f> < 0.01 and 0.12 < <f> < 0.06, re- 
spectively. Two of our transits, 2008-08-29 and 
2008-09-21, coincide with the transit epochs pub- 
lished by S09. The observing log is summarized in 
Tabled! 

The initial trimming, bias and flatficld correc- 
tions of all the collected data were performed us- 
ing custom-made pipelines specifically developed 
for each instrument. The times at the start of 
the exposure are recorded in the image headers, 



2 In the remaining of the text we refer to each individual 
transit by the UT date of mid-time of the transit, using 
the following notation YYYY-MM-DD 



in particular we used the value of the Modified Ju- 
lian Day (JD-2400000.5) field. In the SMARTS 
telescope, the time stamp recorded in the header 
of each frame is generated by a IRIG-B GPS time 
synchronization protocol connected to the com- 
puters that control the instrument. The SOAR 
telescope data use the time values provided by a 
time service connected to the instrument. We con- 
firmed that these values have ~ 1 second precision. 
The time value assigned to each frame corresponds 
to the Julian Day at the start of the exposure plus 
1/2 of the integration time of each image (see sec- 
tion 0] for details). 

WASP-5 is located in a relative empty field, 
where both the target and several well suited com- 
parison stars appear well isolated in our images. 
Therefore, we extracted the flux from the target 
and comparison stars via standard aperture pho- 
tometry, and using our own python-based code. 
We used a range of stellar apertures between 8 and 
12 pixels, and sky rings which extended between 
25 and 35 pixels in radius. 

For each sky-aperture combination, we gener- 
ated differential light curves between the target 
and each comparison star to 1) optimize the aper- 
tures and 2) select the best comparison stars. The 
criterium used in both cases was RMS minimiza- 
tion for the out-of-transit and in-transit data (ex- 
cluding the ingress and egress portions of the light 
curves). The final light curves were generated 
computing the ratio between target's flux and the 
best 2 to 5 comparison stars. 

Finally, some systematics remaining after this 
step were removed by means of linear or quadratic 
regression fits to the out-of-transit light curve 
points using X-Y pixel position, time and/or air- 
mass as free parameters. The final light curves 
present average photometric dispersions of the or- 
der of 0.2% - 0.45%. 

3. Light Curve Modeling 
3.1. Algorithm Comparison 

We performed a comparison between algo- 
rithms that use different statistical uncertainty 
estimation techniques to the transit's parame- 
ters, in order to test potential systematics be- 
tween them. There are different approaches to do 
that statistical error estimation analysis; for exam- 
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pie, JKTEBOFH (|Southworth et alJ l2004albh uses 
the Levenberg-Marquardt Monte Ca rlo (LMMC) 
technique to compute e rrors (see e.g. Southworthl 
2010t Hover et "all 2011 ). while several other stud- 



ies have started to implement Monte Carlo Markov 



Chains (MCMC) techn iques (e.g. lAdams et al 
20ld : iFulton et al.ll201ll ). 



In iHover et al.l (j 2 1 lh we proposed that the re- 
sults of both, the LMMC and MCMC algorithms 
are equivalent if the parameter space lacks of lo- 
cal minima, where LMMC minimization can be 
trapped. Here we further test that proposal by 
comparing the results of both algorithms on the 
WASP-5b data used for this study. We compare 
the results of fitting a light curve of WASP-5b 
with JKTE BOP and the Tra nsit Analysis Pack- 
ag<3 (TAP: lGazak et aL 2011 ). which implements 
the MCMC me thod for the estim ation of errors 
(more details in lFulton et alJlioill and references 
therein). 

Among the parameters that JKTEBOP allows 
to fit are: the planet-to-star radius ratio (R p /R s ), 
the inclination (i) and eccentricity (e) of the or- 
bit, the out-of-transit baseline flux (F ao t), the mid- 
time of transit (T c ), the quadratic limb darkening 
coefficients (//i and ^2), and the sum of the frac- 
tional radii, R = R p /a + R s /a, where R p and R s 
are the absolute stellar and planetary radii, and 
a is the orbital semi-major axis. TAP allows to 
fit for all those parameters except for the latter, 
which is replaced by a/R s . 

For the comparison we left free all the men- 
tioned parameters except a/Rs and R in TAP 
and JKTEBOP, respectively, since they otherwise 
presented convergence problems. We also fixed 
Foot = 1, e = and ^2(1) = and the orbital 
period to P = 1.62843142 days from Fll since 
any variation in this parameter will be detected 
later in our timing analysis. We used 10 4 itera- 
tions in JKTEBOP and 10 chains of 10 5 steps each 
in TAP. We discarded the first 10% iterations on 
each chain to compute the final parameter's values 
and its respective errors. The results on each fit, 
shown in Table [21 reveal that the resultant fit val- 
ues of all parameters common to the JKTEBOP 
and TAP algorithms agree within the error, except 
for 
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Figure [2] shows the distribution of each param- 
eter obtained using the LMMC and the MCMC 
techniques with data from 2008-08-21 transit (sim- 
ilar analysis was done with the other 6 complete 
light curves). From the three bottom panels in the 
Figure[2]it is evident that the la errors (defined as 
the 68% of a Gaussian fit to the parameter value 
distributions) obtained with LMMC are generally 
smaller than those obtained using MCMC, since 
the latter does a more exhaustive exploration of 
the parameter space and therefore performs bet- 
ter error estimations. Also, from the top-panel 
of Figure [5J it can be seen that the LMMC re- 
sults for certain parameters can appear biased 
towards their initial input values. That is the 
case for the linear limb-darkening coefficient, for 
which the value resulting from the LMMC anal- 
ysis is Hi(T) = 0.22 ± 0.12 (the initial value was 
0.296). On the other hand, the distribution of val- 
ues for this parameter on a single epoch as given by 
MCMC does not appear Gaussian, revealing that 
the quality of a single transit in the current data 
does not allow to constrain the values of /ii(/). 
Notice, however, that a Gaussian distribution is 
obtained when fitting several transits simultane- 
ously (see Figure 0] and section 13.21) . 

From the test results above we conclude that 
the LMMC and MCMC techniques arrive to simi- 
lar parameter results. However, because the ap- 
parent underestimation of the errors estimated 
by LMMC we have opted for using TAP for our 
analysis of the full WASP-5b transit dataset and 
the re-analysis of all the available data (see next 
section). This underestimation is due to lack of 
multi-parameter uncertainty estimator and fail- 
ure to account for red noise in the minimization 
(jCarter fc Winnl l2009h as TAP docs. Other ad- 
vantages of TAP include that the code can fit a 
greater number of parameters like linear systcm- 
atics in the datasets, and it allows a simultaneous 
fitting of multiple transits. 

3.2. Final Modeling 

We used TAP to fit the nine new transit light 
curves presented in this paper and all the available 
light curves of the system (seven of Fll, two of 
Dll, two of S09 and the two of A08). 

First, we attempted to model each of the new 
light curves independently, but ran into several 
problems. TAP had difficulties fitting the incom- 
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pletc light curves. Also, when fitting individual 
light curves, parameters such as [i\ did not clearly 
converge to a single value, as already mentioned in 
section I3TT1 and illustrated in Figure [3J To avoid 
these problems we fit the seven complete light 
curves simultaneously, leaving as free parameters 
/Lti(I), i, Rp/R s , T c , Foot, in addition to possible 
linear trends to the light curves, F s i ope , and white 
(uncorrelated) and red (correlated) noise compo- 
nents, <j w and oy, respectively. The orbital period, 
the eccentricity, and the longitude of the perias- 
tron were fixed to the values P= 1.62843142 days 
(the value obtained by Fll), e = and lo = 0. 

We used a quadratic limb-darkening law, but 
found that the precision of the light curves was 
not enough to reliably fit the quadratic coefficient, 
so that value was also fixed to = 0. 32, based 

on the tabulated results in IClaretl ( 2000t ) . 

As mentioned in Section [2j we initially cor- 
rected for systematic trends in the light curves 
using linear or quadratic regression fits. Altough 
slopes in the light curves are not clearly appar- 
ent, we leave Foot and F s i ope as free parameters 
to ensure that any small residuals arc properly fit. 
This might create concerns about wheter this two- 
step fitting of systematics can affect the results of 
the fits. To ensure we are not introducing any 
bias on the determination of the planetary param- 
eters, we fit the two sets of data (i.e. the light 
curves with and without systematics trends re- 
moved) with TAP and arrive to consistent values 
of all derived planetary parameters. 

We also searched for potential parameter cor- 
relations in the light curves using the fit results 
of the 2008-08-29 transit described in the previous 
section, where all the parameters were let to vary. 
The resultant parameter correlations are shown in 
Figure |3l This figure reveals a strong correlation 
between a/R s and i. There is also evidence of 
weaker correlations between those two parameters 
and R p /R s . Therefore, to minimize the impact 
of those correlations in our results, we fixed a/R s 
in all the light curves to 5.37 (from Fll), while 
closely monitoring R p / R s and i for variations. 

To fit the transits we ran 10 MCMC chains 
of 10 5 links each, discarding the first 10% results 
from each chain to avoid bias toward the initial in- 
put values of each fitted parameter. Because the 
resulting MCMC distributions for ^ti(i) are not 
Gaussian (see Figures [2] and [3]), that parameter 



was fit simultaneously for all seven light curves, 
while for the other parameters we obtained one 
value per curve and combine them afterward via a 
weighted average. The resulting average values for 
each parameter are listed in Tabled together with 
their la errors. As an example, Figure [4] shows the 
resultant MCMC distributions of i, R p /R s , and T c 
for the transit observed on 2008-08-29, while the 
distribution of fii(I) correspond to the results of 
the simultaneous seven transits fit. This distribu- 
tion is now clearly Gaussian in contrast with the 
previously obtained. 

Finally, we adopted the values of all the param- 
eters that define the shape of the transit derived 
in the fit above and used them as fixed values in 
the two incomplete light curves (2008-10-22 and 
2009-08-06 transits) to derive their mid-times of 
transit, T c . The Foot, F s i ope , a w and oy are still 
left variable in this case. 

Fll used a procedure based on \ 2 minimization 
for modeling their light curves. We re-analyzed 
their data to do an homogeneous study of all the 
light curves, given that a multi-parameter mini- 
mization based on MCMC is statistically more ro- 
bust. We modeled the seven light curves of F110, 
the two light curves of Dll (data provided by 
the author, private communication), the two light 
curves of SOeH] and the two of A08 (data provided 
by the author, private communication) in a similar 
manner to our complete light curves above. The 
Fll transits were observed with a Bessel I filter, 
the Dll and the S09 with a R filter, and the A08 
with R and SDSS i' filter; therefore, we fit one 
simultaneously for all Fll curves, one fii(R) 
for the Dll curves and one for the S09 curves, and 
separate /J,i(i) and /Ji(i?) for the A08 curves. We 
fixed [i2 = 0.32 in all cases. 

The obtained parameters are summarized Table 
|3l The resultant models to all 22 light curves are 
illustrated in Figure [U 

We point out that the errors of the Fll's light 
curves estimated by us are, in average, 70 % larger 
than the reported by Fll. We checked that the 
origin of this difference was not due only by the 
different red-noise estimator methods. Using the 



5 The data is available in the on-line material from the Fll 

publication on PASJ 
''The data is available at the CDS 

l |http://cdsweb.u^s trasbg.fr/ ) 
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same red- noise factor estimated by Fll, we have 
obtained errors consistent with t hose we present in 
Table ICarter fc Wkml (|2009l ) found that time 
averaging and residual permutation methods un- 
derestimated the errors by 15 — 30% compared 
with the wavelet-based method (implemented by 
TAP). 

Using the model results is possible to look for 
variations in the most relevant parameters, in par- 
ticular i and R v j R s , that can reveal the presence 
of an additional body in the system. In Figure [5j 
we plot R p /R s and i as a function of the transit 
epoch, based in the results of the twenty transit 
fits (our two incomplete light curves were not in- 
cluded). We do not see any significant variations 
in those parameters. The weighted average values 
of i and R p /R s based on all the light curves results 
are summarized in Table |4] We studied in detail 
the timing of the transits in the next section. 

4. Timing Analysis 

The times in our nine transit data and the 
Dll data were initially computed in Coordinated 
Universal Time (UTC) and then converted to 
Baryccntric Julian Days, e xpressed in Ter r estria l 
Time, BJD(TT), using the lEastman et all (|201Ct ) 
online calculator Q. The transit times of S09 and 
A08, which were initially expressed in HJD(UT) 
have also been converted to BJD(TT). No conver- 
sion was applied to the light curves reported by 
Fll. 

The times of the common transits, 2008-08-29 
and 2008-09-21, derived from our light curves are 
consistent within the errors in the values derived 
by us and also by Fll from S09 data. 

Using the Fll's ephemeris equation, we cal- 
culated the residuals of the mid-times of the 22 
transits of WASP-5b analyzed in this work. The 
top panel in Figure |6l shows the Observed minus 
Calculated (O — C) diagram for our nine transits. 
In the middle panel of the figure we combine the 
O — C values of our nine transits with the new val- 
ues derived for the Fll, Dll, S09 and A08 (shown 
as open circles). As illustrated in that figure, a 
linear trend with a slope of 2.54 x 10~ 6 days is ob- 
served in the time residuals of all transits. That 



7 http: / /astroutils.astronomy.ohio- 
state.edu/time/utc2bjd.html 



trend can be explained by the accumulation of er- 
rors in the current orbital period and To of the 
transits over time, and therefore can be modeled 
out. 

This linear regression of the points in the O — C 
diagram has a xled = 1 - 22 (x 2 = 24.37 for 20 de- 
grees of freedom), which is significantly smaller 
than the value obtained for Fll of x 2 red = 3.66 
(x 2 = 32.2 for 9 degrees of freedom). Addition- 
ally, we confirmed that with our results for the 11 
epochs included in Fll's analysis we also obtained 
an smaller x 2 (x 2 = 15.45 that yields xled = 1-72). 
This result lies in the fact that our T c uncertainties 
are larger than those estimated by Fll. 

Once the linear trend is removed the updated 
ephemeris equation is: 



T c = 2454375.62459(23) [BJDtt] + 
1.62842888(78) x E, 

where T c is the central time of a transit in the 
epoch E since the reference time T . The errors of 
the last digits are shown in parenthesis. The bot- 
tom panel in Figure [6] shows the resulting O — C 
values of all available transits using the updated 
ephemeris equation. The resultant O — C dia- 
gram is consistent with a constant period, and we 
conclude that the observed TTV residuals (with 
a RMS of ~ 0.00073 days ~ 63 seconds), are 
most likely introduced by data uncertainties and 
systcmatics rather than due by gravitational per- 
turbations of an orbital companion. This newly 
obtained precision permits to place strong con- 
straints in the mass of an hypothetical companion, 
particularly in MMR's, as we discuss in the next 
section. 

5. Limits to additional planets 

To place upper limits to the potential per- 
turbers in the WASP-5 system based in the de- 
rived TTV RMS of about 60 sec we use Mercury 
( Chambers 19991 ) N-body simulator. The input 
parameters to Mercury include the mass and the 
radius of both the star and the transiting planets, 
the planet-to-star orbital separation, as well as the 
inclination, eccentricity and periastron longitude 
of the system. The values for all these parame- 
ters were adopted from S09. In addition, all the 
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initial relative angles between the perturbcr and 
WASP-5b were set to zero. 

We explored a wide range of perturber masses 
between 1 M ffi and 4000 M® in initial steps of 
50 Mq, which are subsequently refined as de- 
scribed below. For the semi-major axis distances 
we explore a range between 0.001 and 1.2 AU in 
steps of 0.001 AU, which was further reduce near 
resonances. The density of the perturber was kept 
constant to that of Earth for M p < 10M® and to 
that of Jupiter for Mp > 200 M®, it was varied 
linearly for masses in between. Also, we assumed 
the perturber to be in a circular orbit and coplanar 
to WASP-5b, since this configuration provides the 
most strict limit to the amplitude of the TTVs 
for a given perturber's mass. Non-zero eccen- 
tricities and non-coplanar orbits p roduce large r 
TTVs as alr e ady p ointe d out by e.g. Bear] (120091) , 



Hover et al.l (|201lh and iFukui et al.l (|201lh . For 



each model configuration we let the system relax 
for five years, and then we used the next five years 
to obtain our fit results, which in total is more 
than 3 times the time span of the observations. 
These 5 years of relaxation time permits to mini- 
mize the effect of any initial bias (e.g. the relative 
angles). We found orbits between 0.02 and 0.035 
AU to be unstable due to the presence of WASP- 
5b. For all other (stable) orbits we recorded the 
central times of each transit of WASP-5b and com- 
puted the predicted TTVs for each configuration, 
assuming an average constant period. Addition- 
ally, we checked that the fitted average period did 
not deviate by more than 3 a from the obtained 
orbital period of WASP-5b. Also, to ensure a good 
sampling of the potential perturber's mass, we re- 
duced the steps in M pert to IMq whenever the 
TTVs approached 60 sec. 

The results of our model simulations is illus- 
trated in Figure [3 where we show the M pert (M®) 
versus a(AU) diagram that places the mass limits 
to potential perturbcrs in the system. The solid 
line in the diagram indicates the derived upper 
limits to the mass of the pcrturbers that would 
produce TTVs RMS of 60 sec at different orbital 
separation. The dashed line shows the perturbcr 
mass upper limits imposed by the most recent ra- 
dial velocity observations of the WASP-5 system, 
for which we have adopted a precision of 15 m/s 
(A08 and lTriaud et al-lboiol report RV precision 



Figure [7] thus shows that the perturber would 
have been detected by RV measurements in all 
areas except around the 1:2, 5:3 and 2:1 MMRs, 
where it could have a maximum mass of 5, 1 and 
2 M®, respectively. 

6. Conclusions 

We present nine new transit light curves of 
WASP-5b. We homogeneously model these light 
curves together with all available transit data of 
this system. Based in these fits we search for any 
variation in the timing of the transits. 

Using 22 transit epochs we updated the ephemeris 
equation and we find a TTVs RMS of 63 seconds. 
All the transit times are consistent with a constant 
orbital period within 2a. 

Our linear fit of the transit times has a 
Xreduce = 1-22, whic h is considerably lo wer than 
the value found by IFukui et"aL ( 201lh used to 
implied the presence of an perturber body. 

Despite obtain ing a similar TTV RMS than 
Fukui et al.1 (|201ll ) (~ 1 min), we conclude a much 
smaller significance to deviations from a constant 
period due to our larger per-epoch uncertainties 
as obtained by the MCMC algorithm. 

If the system has an additional orbiting body, 
its mass has to be lower than 5, 1 and 2 M®, in the 
2:1, 5:3 and 1:2 resonances. In any other location 
the perturber would have been detected by RVs. 

We search for any trend in the depth of the 
transit and inclination of the orbit but we do not 
see any clear evidence of variation with statistical 
significance. 
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of 14 m/s and 12 — 18 m/s, respectively). 
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Table 1 

Observational information of each night. 



Transit Date 


Telescope /Instrument 


Filter 


Integration Time [s] 


airmass range 


Epoch 


2008-08-21 


SMARTS-lm/Y4KCam 


Bcsscll 1 


13 


1.7 - 


1.01 


199 


2008-08-29" 


SMARTS-lm/Y4KCam 


Bessell I 


10 


1.05 


- 1.02 - 1.06 


204 


2008-09-21" 


SMARTS-lm/Y4KCam 


Bcsscll 1 


10,7 


1.9 - 


1.01 - 1.07 


218 


2008-10-22 b 


SOAR/SOI 


Bcsscll 1 


7,5,3 


1.12 


- 1.02 


237 


2008-11-04 


SOAR/SOI 


Bcsscll 1 


3 


1.07 


- 1.02 - 1.4 


215 


2008-11-17 


SMARTS-lm/Y4KCam 


Bcsscll I 


10 


1.02 


- 1.4 


25.3 


2009-06-22 


SOAR/SOI 


Bcsscll I 


7,5,3 


1.95 


- 1.02 


387 


2009-08-06 b 


SOAR/SOI 


Bcsscll 1 


5,4 


1.07 


- 1.15 


414 


2009-10-25 


SMARTS-lm/Y4KCam 


Bcsscll I 


15 


1.06 


- 1.02 - 1.97 


463 



a This transit was also observed byfSoulhworth ct al. (2009). 
b This transit has a incomplete phase coverage. 
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Table 2 

Values obtained with Levenberg-Marquardt Monte Carlo (JKTEBOP) and Markov Chain 
Monte Carlo (TAP) algorithms with data of the 2008-08-21 transit of WASP-5b. 



Parameter 


JKTEBOP 


TAP 


Rp/Rs 


0.0988 ±0.0018 


0.0988 ±0.0026 




83.4 ± 1.5 


83.7 ±2.3 




0.22 ±0.12 


0.45 ±0.11 


T c - 2454699 (UT) 


0.67690 ± .00035 


0.67697 ±0.00041 


(R p + R s )/a 


0.223 ±0.015 




a/R a 




5.01 ±0.48 
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Table 3 

Adjusted parameters for each transit using TAP. 



Transit date Epoch 



R P /R S 



(BJDtt) 




+0.00040 


1.6 
2.8 


'-0.00041 
+ 0.00052 


-0.00051 
+ 0.00075 


2.!) 


-0.00076 
+ 0.00047 


2.:! 
5.3 
2.5 
10.3 
8.2 


'-0.00045 
+ 0.00030 


'-0. 00030 
,+0.00069 


-0.00066 
+ 0.00036 
-0.00036 
+0.00080 


-0-00079 
+ 0.00042 


5.6 


-0.00043 



2008-08-21 


199 





2008-08-29 


204 





2008-09-21 


218 





2008-10-22 


2:S7 




2008-11-04 


245 





2008-11-17 


25.3 





2009-06-22 


387 





2009-08-06 


414 




2009-10-25 


463 






I i o+O.OOlS 

II -0.0015 

i n9+ 0019 

U -0.0020 
QgQ+0.0027 
UOU -0.0026 

0.1116 b 



I , r+0.0027 
LJ - -0.0028 
I n1 +0.0022 
U -0.0024 

0.1116 b 

+ 0.0020 
-0.0021 



0.1114 



85.60±S:al 

o= r q +0.2S 

85 7fi+°- 46 
00 - '"-0.40 

85.47 b 

+ 0.17 
0.16 
+ 0-36 
0.33 
0-21 
0.20 

85.47 b 

+ 0.26 
-0.23 



85.17 
85.45 
85.62 



5.71! 



0.237 
0.237 
0.237 


0.237 
0.237 
0.237 


0.237 



+ 0.05 

0.049 
+ 0.05 

0.049 
+ 0.05 
0.049 

24 b 

+ 0.05 

0.049 
+ 0.05 

0.049 
+ 0.05 
0.049 

24 b 

+ 0.05 
0.049 



4699.1 
4707.1 
4730.1 
4761.! 
4774.! 
4787.1 
5005.! 
5049.' 
5129.! 



2008-06-18° 


160 


n , , „, +0.0032 
U -- L1 ' i -0.0032 


85.02l° 44 


0.2921°, 


089 
089 


4636.174591° 


00079 
00082 


2.!) 


2008-11-02° 


244 


1109+ 0034 

u. 11U3_q 0032 


85.62+°/™ 


292+"" so 


4772.962121° 


(10071 
.00075 


2.2 


2009-09-04° 


432 


1095+ 0048 

U. -LU»U_ 0047 


85.54j:°;« 


0.292+° 


080 
(ISO 


5079.108301° 


00075 
.000711 


2.0 


2009-10-05° 


451 


1091+ 0041 

u. 1U31_„ 004B 


85.44±S;»S 


0.2921° 


080 
089 


5110.046071° 


0008 7 
.000811 


10.1 


2009-10-18° 


459 


1096+°°° 3 ° 

U - 1UJD -0.0031 


86.13«:f 7 


0.2921° 


089 
089 


5123.076111° 


1)0070 
0007!) 


2.5 


2010-06-16° 


607 


112 i+ 0044 

U. J-1_1_q QQ42 


87.30+q| 8 


0.2921° 


080 
089 


5364.08151°, 


0011 
0011 


4.9 


2010-06-29° 


615 


1097+ 0040 

U.1UJ(_ 00 44 


85.67lS; 6 4 3 8 


0.2921° 


080 
080 


5377.109551°, 


00001 
.00093 


5.2 



2009- 09-01 d 

2010- 09-09 d 



430 
659 



111lt"-WH» 
-0 .0029 

ii54+°-°° 41 



86.16 
85.92 



0.59 
0.53 
+ 0.94 
0.68 



0.51 
0.51 



+0.11 

0.13 
+0.11 
0.13 



5075.849471™ 

7 + 0.0010 
-0.0011 



5448.75927+ 



8.0 
4.6 



2008-08-29° 
2008-09-21° 



204 
218 



1109+ 0011 

u - llu3 -0 .0010 

1102+°°° 14 
u.j-j.u - .o 015 



85.78l° ,f 8 
85.78l°; 2 2 4 



+ 0.052 
0.053 



0.367 

U - JO '_0.053 



Avn-7 oocoQ + 0-00023 

4707.82523_ 00025 

4IOU.D^M_ .00031 



3.0 

3.7 



2007-10-10 1 
2007-10-13*' 



0.1095! 



- + 0.0017 
-0.0020 

noi+° °° 61 

u '- Llu - L -0.0066 



85.61 
84.95 



+ 0.37 
0.29 
0.5 9 



0.37Io i 
39+°- 18 



4383.767501™ 
4387.0. 22751°°^° 



4.!) 
11.7 



a Iii all the fits the quadratic coefficient was fixed to {12 — 0.32. 

b These parameters were fixed in the modeling and correspond to the weighted average of the results of the other 
seven full phase covered light curves presented in this work. 

c,d,e,/ Fitti resu it s f the transits of IFukui et al.l II2011D , IDragomir et al.l (I2011fl . ISouthworth et al.1 J2009D and 

Anderson ct al. (2008), respectively. 
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Table 4 

Improved orbital values derived from the weighted average of the light curve's fits. 



Parameter Adopted Value la Error 



a/R s a 
Rpl R s 

Period [days] 
T [BJD TT ] 2454375.62549 ±0.00023 



5.37 
0.1111 
85.56 
1.62842888 



±0.15 
±0.0005 
±0.07 
±0.00000078 



a Va lue and error adopted from iFukui et al 
( 20111) . 



12 



1.0 - 



0.9 ■ 



0.8 



8 



+ 0.7 ■ 

H 

e; 



0.6 



2008-08-21 
2008-08-29 
2008-09-21 
2008-10-22 
2008-11-04 

2008- 11-17 

2009- 06-22 
2009-08-06 
2009-10-25 

2008-06-18 

2008- 11-02 

2009- 09-04 
2009-10-05 

2009- 10-18 

2010- 06-16 
2010-06-29 



This work 





Fll 



0.5 : 



2009- 09-01 

2010- 09-09 



Dll 



2008-08-29 
2008-09-21 



S09 



0.4 



2007-10-10 
2007-10-13 




AO 8 



-0.2 



-0.1 



0.0 



0.1 



0.2 



Fig. 1. — Light curves of the nine transits of WASP-5b presented in this work, the seven transits of Fll, the 
two transits of Dll, S09 and A08. The solid lines show our best model fits using TAP (see section |3~2]) . The 
UT date is indicated in the left of each light curve. 
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,1 
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-0.2 



0.0 



0.2 



0.4 



0.6 



0.8 



Mi) 



80 



82 



84 



86 



88 



90 



i (deg) 



800 



400 



400 " 1 1 1 1 f1 iT'i^Ntom 



0.093 0.095 0.097 0.099 0.101 0.103 

R v /R s 



-^wt-iTTI.I II II II 11 II II r^r*» 



0.6760 



0.6765 



0.6770 

T c -2454699 days (UT) 



0.6775 



0.6780 



Fig. 2. — Histograms of the 10 000 LMMC iterations with JKTEBOP (gray histograms) and of the 10 chains 
with 10 000 links (after discarding the first 10% results on each chain) obtained with TAP (black histograms) 
for the 2008-08-21 epoch. For comparison, the binning factor of the TAP results histograms is 9 times the 
binning factor used for the JKTEBOP results. In both fittings, all the parameters were left free except for 
quadratic limb darkenning coefficients (fi2(I) = 0), eccentricity (e = 0) and the pcriastron longitude (lu = 0). 
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Fig. 3. — Results of the Markov Chain Monte Carlo iterations resulting of fitting the 2008-08-29 transit data 
with TAP, which show the correlation between the fitted light curve parameters a/R s , R p /R s , T c , i and 
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0.11 



0.16 



0.21 



0.26 



0.31 



0.36 




0.104 0.106 



0.108 0.110 

R v /R a 



0.112 0.114 



0.8177 0.8182 



0.8187 0.8192 0.8197 

T c -2454707 days (UT) 



0.116 




0.8202 0.8207 



Fig. 4. — Histograms of the Markov Chain Monte Carlo iterations resulting of fitting the 2008-08-29 transit 
data with TAP, for their fitted parameters: ^i(I), R p /R s , i and T c . The dashed lines show the fitted value 
and the ±lcr errors. 
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Fig. 5. — Derived values of the orbital inclination (Top Panel) and planet-to-star radii ratio, R p /R s 
(Bottom Panel) for all modeled transits. The solid circles correspond to our seven complete transits. The 
open circles correspond to the seven transits of Fll and the two transits of Dll, S09 and A08. The weighted 
average to all points is represented by the solid line on each panel and the dashed lines show the Ha errors 
of those fits. No significant variations are apparent for these parameters in the time span of the observations. 
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Fig. 6. — Top-panel: Observed minus Calculated diagram of the central times of the nine transits reported 
in this work. Middle-panel: O — C residuals of our nine transits (solid circles) combined with the O — C 
residuals of the new fits to the Fll, Dll, S09 and A08 transits (open circles). The dashed line shows the 
linear trend found in the data. Bottom-panel: O — C residuals of all available data after removing the 
linear trend. The solid and point-dashed lines in this figure correspond, respectively, to our new ephemeris 
equation fit and its associated ±lcr errors. 
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Fig. 7. — Upper mass limits of an orbital perturber derived by dynamical simulations done with Mercury 
code (jChamberdll999l ). The solid line represents transit timing variations limits with a RMS of 1 minute. 
The dashed line corresponds to the limits imposed by the current radial velocities observations. Vertical 
lines indicate the location of the MMR distances with WASP-5b for orbital separations of less than 0.08 AU. 
The gray band indicates the range of distances in which any other object would be in a unstable orbit. 



19 



